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I  INTRODUCTION 


This  document  presents  the  final  results  of  an  effort  to  develop  a 
computer  model  of  the  transionospheric  radio  communication  channel.  It 
contains  an  overview  of  the  model,  including  the  propagation  theory  em¬ 
ployed,  the  ionospheric  morphology  incorporated,  and  the  signal  statistics 

★ 

characterized.  Its  companion  volume  (Rino  et  al. ,  1977)  will  be  a  users' 
guide  to  the  computer  codes  developed  in  the  program. 

Initial  work  on  the  program  extended  a  preexisting  model  of  inten¬ 
sity  scintillation  to  provide  a  mathematical  framework  for  characterizing 
the  full  first-order  statistics  of  a  complex  signal  transmitted  through 
the  ionosphere.  In  the  first  half  of  the  present  contract,  the  first- 
order,  signal-statistical  model  was  tested  against  observations  of  inten¬ 
sity  and  phase  scintillation  from  the  Navy  Navigation  (Transit)  satellites 
and  ATS-6.  The  model  was  thereby  improved  substantially  and  extended  to 
the  second-order,  temporal-statistical  domain  to  permit  calculation  of 
the  temporal  autocorrelation  functions  of  intensity  and  phase.  In  addi¬ 
tion,  conditions  were  removed  that  had  restricted  the  model  to  weak, 
single  scatter. 

One  aspect  of  recent  work  on  the  program  was  to  provide  the  capability 
for  calculating  second-order  signal  statistics  to  other  than  the  temporal 
domain.  In  particular,  one  may  now  employ  the  computer  code  IONSCNT  to 
calculate  the  spatial  autocorrelation  functions  of  phase  and  intensity 
and  the  intensity  correlation  coefficient  for  two  frequencies  transmitted 
through  the  ionosphere.  The  final  work  was  to  perform  iterative  testing 
of  the  model  against  definitive  intensity  and  phase  scintillation  data 
obtained  by  means  of  the  DNA  Wideband  satellite,  to  improve  quantitatively 
both  the  signal-statistical  and  morphological  aspects  of  the  model. 


References  are  listed  at  the  end  of  this  report. 


The  signal- statistical  postulates  and  the  propagation  theory  employed 
are  described  in  Section  II  of  this  report.  The  average  morphology  of 
the  ionospheric  structural  parameters  that  control  scintillation,  as  in¬ 
corporated  in  the  model,  is  described  in  Section  III.  The  two  computer 
codes  by  means  of  which  the  model  is  implemented  are  outlined  in  Section 
IV,  and  some  sample  calculations  are  presented;  the  codes  will  be  de¬ 
scribed  fully  in  a  companion  volume  (Rino  et  al.,  1977).  It  is  to  be 
noted  that  the  morphology  employed  in  the  model  is  limited  to  "average" 
conditions  representative  of  the  ionosphere  for  a  set  of  state  parameters 
(location,  time  of  day,  day  of  the  year,  sunspot  number,  and  planetary 
magnetic  index)  to  be  specified  by  the  user.  Departures  from  the  "mean" 
behavior  for  a  given  set  of  state  parameters  are  not  accounted  for.  This 
and  other  limitations  of  the  model  are  described  in  Section  V  of  this 
volume,  in  which  we  present  our  conclusions  from  this  work. 


II  SIGNAL  STATISTICS  AND  PROPAGATION  THEORY 


A.  Heuristic  Description 

The  transionospheric  communication  channel  may  be  regarded  as  a 
linear  filter,  the  transfer  function  of  which  varies  temporally  and  spa¬ 
tially.  Fundamentally,  the  ionosphere  behaves  as  a  dispersive,  phase- 
shifting  network  for  all  radio  waves.  The  gross  carrier  phase  shift, 
group  delay,  and  dispersion  introduced  all  are  proportional  to  the  inte¬ 
gral  of  electron  density  along  the  principal  raypath  of  a  communication 
link.  This  integral,  the  "total  electron  content"  (TEC),  displays  trend- 
like  variations  with  location  on  the  earth,  time  of  day,  season,  and  the 
level  of  solar /geophysical  activity.  These  trends  take  place  on  very 
large  spatial  scales  (hundreds  and  thousands  of  kilometers)  and  long  time 
scales  (tens  of  minutes  and  hours)  ,  and  the  communication-channel  effects 
associated  with  TEC  are  best  treated  deterministically.  Models  of  TEC 
behavior  and  effects  (e.g.,  Klobuchar,  1975)  are  available  for  applica¬ 
tion  to  the  design  and  operation  of  communication  and  navigation  systems 
and  of  radars. 

The  present  work  addresses  effects  of  smaller-scale  (tens  of  meters 
to  tens  of  kilometers)  ionospheric  structure  on  link  performance  over 
shorter  time  scales  (fractions  of  a  second  to  minutes).  Unlike  the  de¬ 
terministic  effects  of  TEC,  the  intensity  and  phase  scintillations  pro¬ 
duced  by  the  small-scale  structure  are  best  described  in  terms  of  signal 
statistics.  As  described  in  interim  reports  on  this  contract  (e.g., 
Fremouw  and  Rino,  1976),  the  signal-statistical  model  developed  is  com¬ 
plementary  to  existing  TEC  models,  in  terms  of  time  scale  and  spatial 
scale  and  in  terms  of  link  effects. 

The  large-scale  structures  identified  as  trends  in  TEC  produce  pure 
phase  trends  in  CW  signals,  whereas  the  unambiguous  signature  of  smaller- 
scale,  scintillation-producing  structure  is  well-developed  intensity 


fluctuation,  accompanied  by  phase  variations.  The  slowest  intensity 
scintillations  typically  are  observed  to  be  relatively  weak  and  to  be 
accompanied  by  relatively  strong  phase  scintillations. 

An  example  of  complex-signal  scintillation  is  presented  in  Figure  1. 
The  data  presented  are  from  the  DNA  Wideband  satellite  (Fremouw  et  al.. 
1977).  Phase  (TEC)  trends  with  Fourier  periods  greater  than  10  s  have 
been  removed,  as  have  intensity  trends  (due  to  changing  range  to  the 
satellite,  etc.).  The  complex- signal  scintillation  clearly  is  dominated 
by  phase  variations,  but  intensity  fluctuations  also  are  present  on  all 
Fourier  frequencies  greater  than  0.1  Hz,  as  is  demonstrated  in  Figures 
2  and  3. 

Figure  2  illustrates  the  scintillations  contained  in  Figure  1  that 
have  Fourier  periods  between  2.5  and  10  s,  isolated  by  means  of  filtering 
intensity  and  phase.  The  weak  intensity  variations  present  are  correlated 
with  the  strong  phase  variations  that  dominate  the  complex-signal  scintil¬ 
lation  in  this  period  range.  The  correlation  is  evident  in  the  complex- 
plane  scatter  plot  at  the  top  of  the  figure;  that  is,  there  is  a  marked 
tendency  for  a  phase  advance  to  be  accompanied  by  a  weak  fade  and  for  a 
phase  retardation  to  be  accompanied  by  weak  signal  enhancement.  This  is 
the  behavior  that  would  be  expected  in  the  geometric-optics  (near)  zone 
of  large,  lenslike  ionospheric  irregularities. 

The  scintillations  with  Fourier  periods  shorter  than  2.5  s  are  illus¬ 
trated  in  Figure  3,  which  shows  that  most  of  the  intensity  scintillation 
occurs  in  this  period  range,  accompanied  by  relatively  weak  phase  fluctu¬ 
ations.  The  statistical  behavior  illustrated  is  consistent  with  diffrac¬ 
tive  scatter  by  smaller-scale  irregularities. 

The  behavior  typified  by  Figures  2  and  3  prompted  us  to  postulate 
that  the  complex  envelope,  E,  of  a  scintillating  signal  may  be  character¬ 
ized  in  terms  of  the  following  two- component ,  signal-statistical  model 
(Fremouw  and  Rino,  1976): 


E  =  Es  Ef  =  (xg  +  iyg)  exp  (Xf  +  i0f) 


(1) 


THREE  REPRESENTATIONS  OF  VHF  (138  MHz)  COMPLEX-SIGNAL 
FLUCTUATIONS  DURING  90  s  OF  A  SATELLITE  PASS  OVER  STANFORD 
CALIFORNIA.  Top  to  bottom:  Scatter  plot  on  complex  plane;  real  and 
imaginary  parts  of  complex  signal;  amplitude  and  phase  of  complex  signal. 
Intensity  and  phase  scintillation  indices  are  S4  =  0.28,  ana  oq  3  1.58  rad, 
respectively. 


TIME  —  s 
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FIGURE  2  THE  FOCUS  COMPONENT  (containing  intensity  and  phase  fluctuations  with 

Fourier  periods  shorter  than  2.5  s)  OF  THE  COMPLEX-SIGNAL  SCINTILLATION 
PRESENTED  IN  FIGURE  1 


It  is 


where  xg  and  yg  are  jointly  Gaussian  variates,  as  are  Xj  and  0^. 

supposed  that  E  and  E  are  statistically  independent.  Thus,  the  first- 
s  I  2 

order  statistics  of  E  are  totally  specified  by  six  (co)variances :  a  , 

2  2  2  x 
a  .  a  ,  a  ,  a.,  a  The  "six  sigmas"  are  related  to  structural  param- 
y»  xy»  0‘ 

eters  of  the  ionosphere  by  means  of  the  propagation  theory  described  in 
Section  II-B;  the  resulting  relationships  are  summarized  in  Section  III-A. 

There  have  long  been  two  intuitive  notions  about  the  complex-signal 
statistics  to  be  expected  for  a  scintillating  signal.  The  first  is  based 
on  the  conceptual  model  of  summation  of  "randomly"  phased  contributions 
from  a  large  number  of  scattering  centers  (Ratcliffe,  1956).  This  notion 
leads  to  the  suggestion  of  Rice  (1945)  statistics  or  a  generalization 
thereof  (Beckman  and  Spizzichino,  1963),  in  which  the  quadrature  components 
of  the  signal  are  possibly  correlated  Gaussian  variates  with  possibly 
different  variances  (Rino  and  Fremouw,  1973;  Rino  et  al.,  1976).  The 
second  notion  is  of  a  radio  wave  propagating  in  essentially  straight  paths 
along  which  the  phase  is  randomly  increased  and  decreased  in  such  a  way 
that  the  output  phase  is  a  Gaussian  variate  (Bramley,  1955).  A  natural 
extension  of  this  view  is  to  suppose  that  the  real  and  imaginary  parts 
of  the  logarithm  of  the  signal  are  normal  variates. 

The  two-component  model  defined  above  is  simply  a  unification  of  the 
foregoing  two  notions  in  a  way  suggested  by  the  first  definitive  complex- 
signal  scintillation  data  to  become  available.  The  two  components  are 
multiplicative  because  they  were  isolated  experimentally  by  a  process 
equivalent  to  a  coherent  AGC  circuit.  Heuristically ,  the  components 
might  be  expected  to  be  multiplicative  because  they  stem  from  additive 
perturbations  of  the  imaginary  part  of  the  complex  logarithm  (i.e., 
phase).  A  partial  formalizing  of  the  two-component  propagation  theory 
is  described  in  Section  II-B. 

Our  propagation  calculations  employ  (1)  the  phase-screen  approach 
to  describing  the  wave  perturbation  imposed  by  ionospheric  structure, 
subject  to  the  Markov  approximation,  and  (2)  the  parabolic-wave  equation 
for  calculating  diffraction  effects,  which  is  appropriate  for  all  narrow- 
angle  scatter  problems.  The  unique  feature  of  the  approach  is  partitioning 


of  the  ionospheric  spatial  spectrum  into  two  scale-size  regimes  (of  no 
geophysical  significance)  to  permit  efficient  computation  and  to  incor¬ 
porate  the  empirically  derived  two-component  signal-statistical  model. 

In  essence,  the  computations  are  divided  into  one  that  accounts  fully  for 
diffractive  scatter  by  relatively  small-scale  irregularities  and  one  that 
employs  the  ray-optics  approximation  appropriate  for  efficient  evaluation 
of  focusing  and  defocusing  by  larger-scale  irregularities. 

A  heuristic  picture  of  the  two-component  model  might  be  of  a  wave 
encountering  a  scattering  layer,  emerging  with  a  scattered  and  a  non- 
scattered  part,  and  then  passing  through  a  region  of  larger  phase-perturbing 
structure.  The  phase  of  both  the  scattered  and  nonscattered  parts  would 
be  spatially  modulated  by  the  second  region.  Statistically,  at  least, 
the  results  would  be  the  same  if  the  positions  of  the  two  perturbation 
processes  occurred  in  the  same  region.  In  this  heuristic  view,  one  might 
suppose  that  only  the  coherent  (nonscattered)  part  of  the  scatter  compo¬ 
nent  would  produce  focuses  after  propagating  through  the  large-scale, 
phase- perturbing  screen.  There  is  some  evidence  to  support  this  suppo¬ 
sition  in  data  from  the  DNA  Wideband  satellite,  and  we  have  imposed  a 
focusing  limit  in  the  model. 


B.  Mathematical  Development 

To  develop  a  mathematical  formalism  that  will  allow  us  to  •  late  the 
parameters  that  characterize  the  first-order  signal  statistics  to  the 
irregularities  that  ultimately  cause  the  scintillation,  we  first  consider 
the  parabolic-wave  equation. 
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By  making  the  partition 


(2) 


an  =  An  +  an, 

s  f 


(3) 


in  Eq.  (2),  where  An  gives  rise  to  the  scatter  component  and  AN  gives 


rise  to  the  focus  component,  and  substituting  E  =  EgE^  from  Eq.  (1)  into 
Eq.  (2),  we  obtain  the  coupled  differential  equations. 
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The  second  equation  is  solved  for  the  log  amplitude,  X^,  and  phase, 
0^,  of  by  using  the  ray  optics  approximations 
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s  Kj, 


The  notation  used  here  is  defined  in  Rino  and  Fremouw  (1977). 


From  Eqs.  (5)  and  (6)  it  follows  that 
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where  R0^  is  the  phase  autocorrelation  function.  The  coefficients  a^ 
are  defined  in  the  companion  volume,  as  is  the  form  of  the  phase  auto¬ 
correlation  function. 

The  scatter  component  is  more  difficult  because  the  near-zone  approxi¬ 
mation  cannot  be  made.  Moreover,  we  must  also  consider  the  coupling  term 
in  Eq.  (4a).  We  assume,  however,  that  the  partitioning  in  Eq.  (3)  is 
such  that  the  scale  for  changes  in  V^E^/E^  (i.e.,  E^)  is  much  larger 

than  the  scale  for  changes  in  7  E  .  Then  Eq.  (4b)  can  be  solved  as  if 

T  S 

the  vTEf/Ef  term  were  constant.  In  this  case,  it  can  be  shown  that  the 

gradient  term  makes  no  contribution  to  the  second-order  moments  of  Eg. 

If  we  make  the  further  assumption  that  Eg  is  Gaussian,  the  second-order 

moments  completely  specify  the  statistics  of  Eg.  If  follows  that  the 

V^E^/E,  term  has  no  effect  on  the  statistics  of  E  --that  is,  E  and  E. 

T  f  f  s  ’  s  f 

are  statistically  independent. 

By  this  argument,  we  can  neglect  the  coupling  term  in  Eq.  (4a).  To 
then  solve  Eq.  (4a),  we  use  the  phase-screen  approximation  so  that 

Es  (P*  z)  =  JJF(P  -  ?')  exp  j  i  0g  (p')  j  dp'  (10) 


where  F(p  -  p  )  is  the  free- space  diffraction  kernel  for  a  wavefield 

scattered  through  narrow  angles.  Because  Eg  is  Gaussian,  its  statistics 

are  completely  specified  by  the  second-order  coherence  functions 

R  *  <E  E  '*>  and  B  ■  <E  E 
s  s  s  s  s  s 

When  R  is  computed  using  Eq.  (10),  the  diffraction  kernel  becomes 

3 

a  delta  function  so  that  R  does  not  vary  with  z.  Thus,  a  simple  analytic 

s 


A 


formula  that  is  given  in  Section  2  of  the  Appendix  of  the  companion  volume 
(Rino  et  al.,  1977),  can  be  used.  For  Bg  a  formula  very  similar  to  Eq. 
(10)  results,  which  must  be  evaluated  numerically.  The  details  are  de¬ 
scribed  in  the  companion  volume. 

To  summarize,  as  long  as  Eg  is  Gaussian  and  the  scale  for  changes 

in  V  in  E  is  very  large  compared  to  the  scale  for  changes  in  V  E  ,  we 
T  f  T  s 

can  relate  the  multiplicative  scatter  and  focus  components  to  disjoint 
and  therefore  independent  spectral  scale  size  regimes  of  the  contributing 
irregularities.  From  these  relations  we  can  calculate  the  signal  moments 
that  describe  the  field  as  a  function  of  parameters  that  characterize  the 
spectrum  of  the  ionospheric  irregularities.  This  constitutes  the  formal 
basis  of  all  the  mathematical  computations  used  in  the  code. 


Ill  SCINTILLATION  MORPHOLOGY 


A.  Controlling  Parameters 

As  will  be  outlined  in  Section  IV,  two  computer  codes  have  been  de¬ 
veloped  under  this  contract.  The  primary  scintillation  calculations  are 
performed  in  the  more  elaborate  of  the  two  codes,  IONSCNT.  There  are 
two  operational  modes  in  IONSCNT,  one  to  provide  first-order  signal- 
statistical  parameters  and  one  to  provide  second-order  quantities.  The 
ultimate  outputs  of  the  first-order  mode  are  the  six  sigmas  that  charac¬ 
terize  the  two-component  scintillation  model  and  the  scintillation  indices 
for  phase,  0,  and  intensity,  I,  which  are  defined  respectively  as  follows: 
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According  to  the  two-component  model,  the  phase  and  intensity  scintil 
lation  indices  for  the  composite  signal  are  given  in  terms  of  those  for 
the  scatter  and  focus  components  as: 
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and 


S*-(S4s  +  l)(S4f  +  1)-1  •  <U> 

The  component  indices  are  obtained  as  follows.  The  most  complicated 
is  the  phase  scintillation  index  for  the  scatter  component.  It  is  ob¬ 
tained  by  performing  the  following  two  integrations: 


over  p  ,  the  Hatfield  phase  distribution  for  a  signal  with  complex  Gaussian 
H 

statistics.  The  Hatfield  distribution  (Hatfield  and  Rino,  1975)  is 
specified  by 


5  [RoS  -  Re  {Boe}] 
I  [Ro„  +  Re  Ks}] 
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where 


R  =  1  -  e 
os 


where  a  is  a  scattering  coefficient  per  unit  depth  in  a  layer  of  thick- 
8 

ness  L  [see  Gq.  (31)],  and  B  is  obtained  from  numerical  evaluation  of 

os 

the  diffraction  integrals  by  a  procedure  described  in  the  companion  vol¬ 
ume  (Rino  et  al.,  1977).  The  quantities  above  are  also  used  to  calculate 
the  intensity  scintillation  index  for  the  scatter  component,  as  follows: 


s?  =  4a2  (i  -  R  )  +  |b  Y 
4s  x  \  os  /  '  os  1 


The  phase  and  intensity  scintillation  indices  for  the  focus  compo¬ 
nent  are  calculated,  respectively,  as 


and 


where  a  L  is  the  variance  of  phase  perturbation  imposed  by  a  layer  of 
t  2 

thickness  L  [see  Eq.  (30)],  and  is  obtained  from  numerical  evaluation 
of  differential  equations  derived  from  the  geometric-optics  approximation 
of  the  parabolic  wave  equation,  as  described  in  Section  II-B;  the  same 
approximation  is  used  to  calculate  aXf0f 

2 

The  propagation  equations  that  govern  ®os»  CTXf»  an<*  CTXf0f  depend  upon 
the  observing  geometry  relative  to  a  horizontally  stratified  ionospheric 
layer  and  to  a  geomagnetic  field  model  and  upon  several  parameters  that 
characterize  the  size  and  shape  of  ionospheric  irregularities.  As  de¬ 
scribed  in  the  companion  volume,  the  three-dimensional  spatial  power  spec¬ 
trum  employed  to  describe  the  irregularities  has  a  power-law  form  with 

a  spectral  index  of  minus  four,  corresponding  to  an  in-situ  (one  dimen- 

_2 

sional)  spectrum  of  K  .  Provision  is  made  for  anisotropy  in  two  di¬ 
mensions,  one  along  the  magnetic  field  and  one  perpendicular  to  the  field. 
The  spectrum  over  which  the  variance  of  ionospheric  electron  density, 

((An)  ),  is  modeled  is  limited  by  an  outer-scale  parameter,  and  the 
spectrum  is  divided  into  two  components  by  means  of  an  inner-scale  cutoff 
for  the  focus  component,  0^,  and  an  outer- scale  parameter  for  the  scatter 
component ,  . 


B.  Behavior  of  Axial  Ratios.  Scale  Sizes,  and  HeiRht 

The  axial  ratios,  a  (along  the  field)  and  b  (perpendicular  to  the 
geomagnetic  meridian  plane)  were  set  to  reproduce  known  geometrical  char¬ 
acteristics  of  scintillation.  Interferometer  observations  have  shown 
that  scintillation-producing  irregularities  are  highly  elongated  along 
the  magnetic  field  near  the  geomagnetic  equator  (e.g.,  Roster  et  al. , 
1966).  Recent  studies  by  Rino  (1977)  indicate  that  high-latitude  irregu¬ 
larities  are  somewhat  elongated  perpendicular  to  the  geomagnetic  meridian 
as  well  as  along  the  magnetic  field  (i.e.,  that  they  are  sheet-like  along 
L-shells  rather  than  axially  symmetric  along  the  field). 
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To  include  the  foregoing  experimental  characteristics,  the  nxial 

ratios  a  and  b  have  been  modeled  as  functions  of  geomagnetic  latitude, 

X  ,  as  follows: 
m* 


and 
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where  X,  demarks  the  geomagnetic  latitude  of  the  scintillation  boundary 

D 

between  the  midlatitude  ionosphere  and  the  subauroral  ionosphere.  Equa¬ 
tion  (24)  provides  an  along-field  axial  ratio  of  30  near  the  geomagnetic 
equator,  switching  to  10  fairly  abruptly  at  a  geomagnetic  latitude  of 
20°.  Equation  (25)  describes  axial  symmetry  at  latitudes  equatorward  of 
the  scintillation  boundary.  Poleward  of  the  boundary  there  is  an  addi¬ 
tional  elongation  in  the  geomagnetic  east-west  direction,  with  a  ratio 
of  6:1  (b  =  6)  relative  to  the  remaining  cross-field  dimension. 

The  three  scale  parameters  O^,  3^,  and  Of g  were  established  by 

iterative  testing  against  signal-statistical  samples  obtained  from  passes 

of  the  DNA  Wideband  satellite.  The  most  satisfactory  fits  were  found  to 

be  different  at  different  latitudes.  (No  data  were  available  from  the 

polar  caps.)  Variation  of  the  scale-size  parameters  with  geomagnetic 

latitude,  X  ,  has  been  modeled  as  follows: 
m 
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Equation  (26)  provides  an  outer  scale  (transverse  to  the  magnetic 
field)  of  just  under  1  km  (i.e.,  approximately  equal  to  the  Fresnel-zone 
radius  for  VHF  waves  in  the  F  layer)  near  the  geomagnetic  equator  and  of 
about  2.6  km  poleward  of  20°  geomagnetic  latitude.  It  is  to  be  noted 
that  the  actual  outer  scale  of  ionospheric  irregularities  probably  is 
considerably  greater  than  one  or  a  few  kilometers;  the  present  model 
accounts  only  for  irregularities  that  produce  scintillation,  defined  in 
terms  of  intensity  fluctuation.  Equations  (27)  and  (28)  are  used  to 
maintain  effective  continuity  between  the  focus  and  scatter  components 
of  the  spatial  spectrum  employed,  while  also  safeguarding  validity  of 
the  ray-optics  approximation  used  for  calculating  propagation  of  the  focus 
component  of  the  perturbed  wave. 

The  reason  for  employing  a  smaller  outer  scale  in  the  equatorial 
region  than  elsewhere  is  that  we  have  consistently  observed  with  the  DNA 
Wideband  satellite  that  the  ratio  of  to  is  larger  (unless,  of  course, 
S4  is  near  its  limiting  value  of  approximately  unity)  in  the  equatorial 
region  than  elsewhere,  for  data  processed  in  an  identical  fashion.  In 
terms  of  the  phase- screen  model,  this  observed  fact  means  that  a  given 
level  of  phase  perturbation  (within  the  spatial  spectrum  admitted  by 
10-s  detrending)  results  in  stronger  propagation  effects  (focusing  and 
diffraction)  near  the  equator  than  at  higher  latitudes.  Employing  a 
smaller  outer  scale  for  a  given  level  of  ((An)  )  is  one  way  of  accounting 
for  such  behavior.  A  systematic  decrease  in  the  spectral  index  with 
decreasing  latitude  also  could  account  for  this  effect.  A  variable  spec¬ 
tral  index,  however,  greatly  complicates  the  computer  code.  Thus,  we 
have  not  incorporated  this  option. 

A  greater  propagation  distance  (i.e.,  employing  a  greater  scattering- 
layer  height)  also  increases  the  level  of  intensity  scintillation.  There 
is  evidence  from  backscatter  observations  (Woodman  and  La  Hoz,  1976),  for 
development  of  fine  structure  from  the  bottomside  to  the  topside  of  the 
equatorial  ionosphere.  We  have  incorporated  a  simple  latitude  dependence 
of  scattering- layer  height  in  the  model,  as  follows: 
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which  places  the  center  of  the  layer  at  350  km  at  geomagnetic  latitudes 
greater  than  about  20°  and  at  500  km  near  the  magnetic  equator. 


C.  Behavior  of  Scattering  Strength 

Once  geometric  and  signal-statistical  behavior  is  modeled  by  choice 
of  axial  ratios,  scale  sizes,  and  height,  the  worldwide  morphological 
behavior  of  scattering  strength  is  the  dominant  consideration  for  scintil¬ 
lation  modeling.  The  total  disturbance  imposed  on  the  wave  consists  of 

the  wavefront  perturbation  imparted  by  both  the  focus  and  scatter  portions 

2 

of  the  spatial  spectrum,  whose  integral  is  ((AN)  ).  For  the  focus  com¬ 
ponent,  the  variance  of  the  phase  perturbation  imposed  by  a  layer  of 
thickness  L  is  given  by 
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The  corresponding  quantity  for  the  scatter  component  is 


G  cos  9 


L  sec  9 


(31) 


where  9  is  the  incidence  angle  on  the  layer,  G  depends  on  details  of  the 
scattering  geometry  through  the  anisotropy  of  the  irregularities,  r^  is 
the  classical  electron  radius,  X  is  the  radio  wavelength,  and  the  K's 
denote  the  modified  Bessel  function. 

All  the  geophysical  parameters  in  Eqs.  (30)  and  (31)  were  discussed 
in  Section  III-B,  except  L  and  ((AN)  )  ■  ((An)  )  +  ((AN)  ).  In  conjunc- 

S  I 

tion  with  the  latitudinal  dependence  of  layer  height,  the  layer  thickness 
also  was  modeled  as  a  simple  function  of  geomagnetic  latitude,  as 
follows: 


>  '  t 
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X  -  20°  \ 
erf  m  - J  km 


With  a  and  ot  set,  the  relative  amount  of  perturbation  assigned  to  the 
I  s 

focus  and  scatter  components  is  that  required  for  their  respective  spectra 
to  be  asymptotically  continuous.  The  condition  is 


-  T  (<“>’)  • 


Finally,  observed  morphological  characteristics  of  intensity  and 

phase  scintillation  were  invoked  to  write  the  following  mean,  morphological 

2 

description  of  ((AN)  ): 


<<■>!)  ■  M„(X.*T)  +  A"h(\.>T-Kp)  +  "‘i(;VT-Kp) 
+  ANe(Xn,Xg,T,",R)  . 


The  first  two  terms  and  the  last  term  on  the  right  of  Eq.  (34)  account 

for  behavior  at,  respectively,  middle,  high,  and  equatorial  latitudes, 

as  functions  of  geomagnetic  latitude,  X  ,  geographic  latitude,  X  ,  local 

m  g 

time  of  day  in  hours,  T,  day  of  the  year,  D,  mean  sunspot  number,  R,  and 
the  planetary  magnetic  activity  index,  K  .  The  third  term  in  Eq.  (34) 
accounts  for  scintillation  associated  with  the  auroral  oval. 

The  midlatitude  term,  centered  at  a  geomagnetic  latitude  of  32.5°, 
contains  only  a  simple  diurnal  variation,  as  follows: 

9/  TTCM  fX  -  32 .5°  J2 

An  =  1.3  X  10  1  +  0.33  cos  exp  — - r— ! —  .  (35) 
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The  high- latitude  term, 

AN.  *  3  X  10 9/  1  +  erf  — t)el/m3,  (36) 

"  \  Xh  / 

dominates  poleward  of  the  scintillation  boundary,  the  geomagnetic  latitude 
and  width  of  which  vary  with  magnetic  activity  and  time  of  day  as  follows: 
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Xh  =  5  +  0.2  Xb  deg 
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Superimposed  on  the  high-latitude  term  is  the  following  auroral-oval 
term: 


where 
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The  most  complicated  term  in  the  morphological  model  is  the  following 
one  describing  behavior  near  the  geomagnetic  equator: 
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where 


TR  " 


6  if  T<  22.5 


Equation  (41)  describes  a  region  of  disturbance  that  peaks  some  10°  on 
each  side  of  the  geomagnetic  equator.  This  somewhat  speculative  morpho¬ 
logical  behavior  has  been  invoked  to  account  for  a  very  persistent  feature 
of  scintillation  observed  at  Ancon,  Peru  and  Kwajalein,  Marshall  Islands 


(both  within  a  few  degrees  of  the  magnetic  equator)  on  signals  from  the 
DNA  Wideband  satellite.  The  persistent  behavior  is  an  increase  in 
scintillation  near  the  ends  of  satellite  passes  that  is  stronger  than  we 
have  been  able  to  account  for  on  the  basis  of  any  purely  geometrical 
feature  of  the  model. 

Equation  (41)  also  includes  a  linear  increase  of  scintillation  activ¬ 
ity  with  increasing  sunspot  number,  and  a  diurnal  peak  somewhat  before 
local  midnight,  with  activity  increasing  rather  rapidly  after  F-layer 
sunset  and  decaying  more  slowly  after  the  peak.  The  seasonal  variation 
includes  equinoctial  peaking  but  with  a  "filling  in"  of  activity  in  the 
local  summer  season  (i.e.,  centered  on  the  December  solstice  in  the  geo¬ 
graphic  southern  hemisphere  and  on  the  June  solstice  in  the  geographic 
northern  hemisphere).  This  somewhat  complicated  behavior  is  consistent 
with  differences  in  seasonal  patterns  observed  at  Kwajalein  and  at  Ancon 
and  reproduces  reasonably  well  the  equinoctial  peaking  observed  by  Foster 
(1966)  and  Foster  et  al.  (1966)  for  observations  at  a  station  near  the 
geographic  equator  in  Africa;  it  is  untested  at  other  locations,  however. 


IV  COMPUTER  CODES 

Two  user  codes  have  been  developed  under  this  contract:  IONSCNT  for 
calculating  first-order  and  second-order  signal-statistical  moments,  in¬ 
cluding  scintillation  indices  for  intensity  and  for  phase  and  several 
correlation  functions;  and  DIST  for  calculating  first-order,  probability- 
density  functions  (PDFS)  for  intensity  and  phase  and  the  corresponding 
cumulative  distribution  functions  (CDFS). 

Figure  4  shows  a  simplified,  overall  flow  diagram  for  IONSCNT.  The 
user  inputs  a  set  of  selected  ionospheric  state  parameters  plus  a  descrip¬ 
tion  of  his  intended  operational  geometry.  The  code  then  calculates  the 
appropriate  ionospheric  structure  parameters  from  the  morphological  model 
contained  in  subroutines  RMSDN  and  MDLPRM.  The  quantities  relevant  to 
the  scatter  and  focus  components  are  passed  to  routines  for  calculating 
diffraction  and  geometric-optics  propagation  effects,  respectively. 

The  outputs  from  IONSCNT  always  include  the  intensity  and  phase 
scintillation  indices,  and  and  the  six  sigmas  that  describe  the 
two- component ,  first-order,  signal-statistical  model.  If  the  user  has 
selected  operating  Mode  II,  second-order  statistical  counterparts  are 
put  out  as  well,  including  the  temporal  autocorrelation  functions  of 
intensity  and  phase,  analogous  information  on  the  spatial  autocorrelation 
function,  and  (if  requested)  the  correlation  coefficient  for  intensity 
fluctuations  on  two  user-selected  frequencies. 

The  temporal  autocorrelation  functions  are  computed  from  the  calcu¬ 
lated  spatial  autocorrelation  function,  using  a  velocity  that  accounts 
both  for  line-of-sight  motion  through  the  F  layer  and  an  assumed  drift 
of  otherwise  unchanging  ionospheric  structure.  The  drift  velocity  is 
eastward  at  a  speed  given  by  the  following  geomagnetic-latitude-dependent 
model: 
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which  yields  50  m/s  in  the  equatorial  zone,  20  m/s  at  middle  latitudes, 
and  upwards  of  100  m/s  poleward  of  the  subauroral  scintillation  boundaries, 
depending  on  the  level  of  geomagnetic  disturbance. 


An  example  of  the  Mode-I  output  of  IONSCNT  is  presented  in  Figure  5. 
The  increase  in  scintillation  early  in  the  second  half  of  the  pass  arises 
from  a  combination  of  (1)  geometric  enhancement  by  sheet-like  irregularit: 
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FIGURE  4  SIMPLIFIED  FLOW  DIAGRAM  FOR  THE  MAIN 
COMPUTER  CODE,  IONSCNT 
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respectively)  appear  in  the  rightmost  two  columns. 


as  the  line  of  sight  grazes  the  L-shell  through  the  sta'.  ior,  and  (2)  .. 

contribution  from  the  auroral  term  in  Eq.  (34).  Graphs  of  S.  and  0 

4  rms 

for  a  representative  pass  of  the  Wideband  satellite  observed  at  Poker  Flat, 
Alaska  under  the  specified  conditions  are  shown  in  Figure  6  for  comparison. 
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FIGURE  6  STANDARD  SUMMARY  DISPLAY  OF  VHF  (138  MHz)  INTENSITY  AND  PHASE 
SCINTILLATION  INDICES  (S4  and  <ty,  respectively)  OBSERVED  IN  ALASKA  ON 
A  PARTICULAR  SATELLITE  PASS  UNDER  CONDITIONS  APPROXIMATED  BY 
THE  INPUT  CONDITIONS  FOR  THE  CALCULATION  ILLUSTRATED  IN 
FIGURE  5 


Note  that  the  pass  illustrated  was  one  of  those  used  in  iterative  testing 
of  the  model;  it  was  chosen  as  one  representative  of  scintillation  in 
Alaska  under  the  specified  conditions  (including  moderate  magnetic  activ¬ 
ity).  A  range  of  scintillation  behavior  presumably  would  be  observed  for 
the  same  set  of  specified  state  parameters,  and  the  model  does  not  account 
for  such  variations. 


In  addition  to  the  scintillation  indices,  the  output  illustrated  in 
Figure  5  gives  relatively  subtle  signal-statistical  information,  in  the 
form  of  the  six  sigmas  (four  variances  and  two  covariances)  that  charac¬ 
terize  the  two  first-order,  signal-statistical  components.  They  appear 
in  the  six  columns  following  the  transmitter  latitude  and  longitude  and 

the  angle  (PSIBP)  between  the  line  of  sight  and  the  geomagnetic  field. 

2  2 

Taking  point  8  in  Figure  5  as  an  example,  we  have  =  0.055,  ay  =  0.271, 
2  2  s 

CTxsys  =  -0.005,  oXf  =  0.000,  a0f  =  6.839,  <?Xf0f  =  -0.002. 

By  comparison,  the  corresponding  values  for  a  30-s  data  segment 
(starting  at  0217:35)  recorded  during  a  geometrically  similar  portion  of 
the  satellite  pass  illustrated  in  Figure  6  are  0.06,  0.27,  -0.02;  0.01, 
4.40,  and  -0.07.  Similarities  between  the  calculated  and  observed  values 
suggest  that  the  existing  model  would  be  useful  for  calculating  the  joint 

pdf  of  intensity  and  phase,  if  the  joint  statistics  should  be  of  interest 

for,  say,  optimal  system  design.  Such  a  capability  has  not  been  incor¬ 
porated  in  the  present  computer  programs  but  could  be  added  to  the  DIST 
code. 

The  DIST  code  now  available  permits  calculation  of  distributions 
for  amplitude  (square  root  of  intensity)  and  phase,  separately.  Examples 

are  presented  in  Figure  7,  using  the  six  sigmas  from  point  8  of  Figure 

5  as  inputs  to  DIST,  As  is  described  fully  in  the  companion  volume  on 
the  codes  (Rino  et  al. ,  1977),  the  method  of  using  DIST  is  to  input  the 
results  of  an  IONSCNT  calculation. 

Figure  8  illustrates  second-order  outputs  from  IONSCNT,  again  using 
point  8  from  Figure  5;  the  format  is  fully  described  in  the  companion 
volume.  The  information  provided  can  be  used  to  characterize  the  spatial 
structure  of  the  field  as  well  as  the  temporal  structure  of  intensity 
and  phase. 

Figures  9  through  14  are  included  to  provide  further  information  on 
the  morphological  behavior  (and  frequency  dependence)  of  scintillation 
as  characterized  by  IONSCNT.  The  reader  will  note  that  the  code  returns 
values  of  S^  slightly  larger  than  unity,  which  is  the  value  corresponding 
to  Rayleigh  statistics,  under  extreme  conditions.  The  two- component , 
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FIGURE  7  EXAMPLE  OF  OUTPUTS  FROM  DIST  CODE,  USING  IONSCNT  RESULTS  FROM  POINT  8  IN  FIGURE  5  AS  INPUTS 


FIGURE  7  (Continued) 
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irregularities  supposed  to  be  aligned  along  L-shells. 
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FIGURE  12  CALCULATED  SEASONAL  VARIATION  OF  NIGHTTIME  SCINTILLATION  AT  ANCON  NEAR  SOLAR  MINIMUM 


lote  differences  near  solstices  compared  with  Figure 


REGION  NEAR  AN  EQUINOX  DURING  A  YEAR  OF  ELEVATED  SOLAR  ACTIVITY.  The  code  would  predict 
Rayleigh  statistics  (S4  =  1.0)  under  these  conditions  for  frequencies  below  50  MHz.  The  phase  scintillation  index 
(PRMS)  varies  approximately  as  f'. 


signal-statistical  model  does  permit  values  greater  than  unity  in  the 
presence  of  substantial  focusing.  The  ray-optics  approximation  employed 
for  calculating  propagation  of  the  focus  component  is  not  always  fully 
satisfactory,  and  a  limit  on  it  has  been  incorporated  in  the  code  for 
conditions  of  strong  scatter.  It  is  to  be  noted  that  values  up  to 
1.2  are  observed  occasionally,  at  least  at  VHF  near  the  geomagnetic  equa 


V  CONCLUSION 


This  report  and  its  companion  volume  (Rino  et  al.,  1977)  describing 
computer  codes  IONSCNT  and  DIST  constitute  the  final  reports  on  this 
contract  to  develop  a  transionospheric  channel  model.  The  work  began  as 
an  attempt  to  collate  previously  existing  and  highly  fragmented  scintil¬ 
lation  data  into  a  form  useful  for  system  application.  During  the  course 
of  the  work,  however,  both  experimental  and  theoretical  work  on  trans¬ 
ionospheric  radio  propagation  progressed  considerably.  The  major  experi¬ 
mental  development  was  the  advent  of  complex-signal  data  from  coherent 
radio  beacons  on  the  Navy  navigation  (Transit)  satellites  ATS-6  and  P76-5 
(the  DNA-002  wideband  beacon).  The  theory  of  propagation  in  a  random 
medium  advanced  both  analytically  and  by  application  of  numerical  calcu¬ 
lation  to  strong-scatter  problems. 

Availability  of  the  coherent-beacon  data  has  yielded  much  insight 
into  the  comp lex- signal  statistics  of  scintillation,  and  this  insight 
has  been  incorporated  in  our  work  to  provide  a  signal-statistical  model 
that  is  complementary  to  deterministic  models  of  TEC  effects.  Our 
ability  to  incorporate  the  new  experimental  findings  has  been  limited, 
however,  by  the  theoretical  framework  we  have  employed.  Signal- 
statistically ,  the  most  fundamental  limitation  of  the  present  model  is 
its  heavy  reliance  on  choice  of  an  outer  scale  to  describe  scintillation- 
producing  ionospheric  structure.  It  is  very  likely  that  the  outer  scales 
invoked  in  the  model  (900  m  in  the  equatorial  region  and  2.6  km  elsewhere) 
have  no  geophysical  significance;  rather  they  represent  the  largest  scales 
that  contribute  to  intensity  scintillation  on  the  paths  employed  in  col¬ 
lection  of  the  data  we  used  for  quantitative  calibration  and  iterative 
testing  of  the  model.  The  paths  employed  were  mainly  between  satellites 
in  high-inclination,  moderate-altitude  (1000-km)  orbits  and  the  ground. 

The  main  caution  to  be  employed  as  a  result  of  the  outer-scale  di¬ 
lemma  arises  in  application  of  the  model  to  phase- sensitive  systems  where 


long-term  stability  is  important,  especially  in  the  equatorial  region. 

The  user  must  always  keep  in  mind  that  the  rms  phase  fluctuations  de¬ 
scribed  by  the  code  are  those  produced  by  ionospheric  structures  on  the 
order  of  a  kilometer  and  smaller  in  east-west  extent  near  the  equator 
(30  km  and  smaller  in  north-south  extent).  Phase  trends  due  to  larger- 
scale  structure  must  be  accounted  for  separately.  For  geostationary 
systems,  this  translates  into  time  scales  substantially  longer  than  20  s. 

Morphologically,  the  most  important  limitation  of  the  model  is  that 
it  does  not  account  for  night-to-night  departures  of  scintillation  con¬ 
ditions  on  a  given  communication  link  during  the  same  season,  at  the  same 
local  time,  and  under  the  same  solar  and  geomagnetic  conditions,  which 
can  be  substantial.  The  model  describes  only  the  "expectation  level"  of 
scintillation  for  a  specified  set  of  observing  conditions.  Data  now 
being  accumulated  could  be  employed  to  include  a  description  of  departures 
from  mean  scintillation  conditions,  but  the  present  model  stops  short  of 
doing  so.  Mean  morphological  limitations  of  the  model  include  some  uncer¬ 
tainty  over  long-term  trends  in  activity  during  a  solar  cycle  and  virtual 
absence  of  any  data  from  the  polar  caps  for  use  in  establishing  the  model. 

Despite  the  foregoing  limitations,  the  model  described  in  this  re¬ 
port  is  the  most  extensive  computer- coded  assimilation  of  scintillation 
information  presently  available  for  systems  planning.  It  will  provide 
systems-relevant  characterization  of  scintillation  for  a  wide  range  of 
operating  circumstances.  Hopefully  it  also  can  serve  as  a  baseline  for 
analysis  of  existing  and  forthcoming  scintillation  data,  in  attempts  to 
understand  the  nature  and  dynamics  of  ionospheric  structure. 
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